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We study the Bose-Einstein condensation for a cloud of 7 Li atoms with attractive nonlocal (finite- 
range) interaction in a harmonic trap. In addition to the low-density metastable branch, that is 
present also in the case of local interaction, a new stable branch appears at higher densities. For a 
large number of atoms, the size of the cloud in the stable high-density branch is independent of the 
trap size and the atoms are in a macroscopic quantum self-trapped configuration. We analyze the 
macroscopic quantum tunneling between the low-density metastable branch and the high-density 
one by using the istanton technique. Moreover we consider the decay rate of the Bose condensate 
due to inelastic two- and three-body collisions. 
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The Bose-Einstein condensation (BEC) for atoms with 
effective attractive interaction, i.e. negative scattering 
length, is an interesting subject of study because of the 
instability of the system: Recent experiments with 7 Li 
vapors in a harmonic trap have shown that the BEC is 
achieved only by a finite and small number of atoms [1]. 
In a uniform system, 7 Li atoms do not form a stable BEC 
state because the s-wave scattering length a s is negative 
and the attractive interaction causes the condensate to 
collapse upon itself [2] . In a harmonic trap, atoms acquire 
zero-point energies which counterbalance the attractive 
interaction: a metastable condensate can be obtained for 
a number of atoms below a critical value N c , which is 
about 10 3 under the condition of present experiments 
[1,3]. 

Recently, we have shown that the inclusion of a non- 
local interaction implies the existence of a new stable 
branch intermediate in density between the known very 
dilute metastable state and the collapsed state [4]. In 
this paper, we discuss the properties of the stable state 
and calculate the tunneling rate of macroscopic quantum 
tunneling (MQT) from the metastable state to the sta- 
ble one by using the istanton technique. Moreover, we 
determine the decay rate of the Bose condensate due to 
inelastic two- and three-body collisions. 

The system of 7 Li atoms is quite dilute and at zero 
temperature the atoms are practically all in the conden- 
sate [1,3]. The action of the system is given by 
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For alkali atoms, the atom-atom interaction is usu- 
ally replaced by an effective zero-range pseudo-potential, 
V(r) — gS 3 (r), where g = 4:nh 2 a s /m is the scattering 
amplitude and a s is the s-wave scattering length. By us- 
ing such an effective potential, and imposing the least 
action principle to Eq. (1), one gets the so-called Gross- 
Pitaevskii (GP) equation [6]. 

In the case of 7 Li, that has negative scattering length, 
the scattering cross section has a significant momentum 
dependence already at very low momenta [7,8]. This mo- 
mentum dependence implies that the effective interaction 
is non local changing from attractive to repulsive at a 
characteristic range r e . We assume [4] that the attrac- 
tive potential has a finite range r e and in addition we 
allow for the presence of a repulsive contribution which 
is modeled as a local positive term defined by a " high en- 
ergy" scattering length a r > 0. The effective interaction 
can be written in the following form: 
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where </>(r, t) is the macroscopic wavefunction of the 
condensed atoms, U(r) is the external potential and 
V{\t — r'|) is the interatomic potential [5]. 



We have considered two choices for the shape function 
f(x): a Lorenzian f(x) = (1 + a; 2 ) -1 and a Gaussian 
f(x) — exp (~x 2 ). The results do not depend on the spe- 
cific choice of f(x) and so we will discuss only the Loren- 
zian case. We use interaction parameters appropriate for 
7 Li (see Parola, Salasnich and Reatto [4]): a s = —27 as, 
r e = 10 3 ob and a r = 6.6 as, where cib is the Bohr ra- 
dius. 

To study the properties of the Bose condensate we 
perform a variational calculation. Because of the ex- 
ternal trap, that we model with an isotropic potential 
U(r) = (l/2)mojQr 2 , we use a Gaussian trial wavefunc- 
tion to describe the Bose condensate and minimize the 
quantum action. In the rest of the paper we write the 
lengths in units of the characteristic length of the trap 
ao = y/ti/ (mujo), the time in units u>q , the action in 
units H and the energy in units Hujq. A good choice for 
the wavefunction is the following 
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where iV is the number of condensed atoms and cr and 
(3 are time-dependent variational parameters, a is the 
width of the condensate. Note that, in order to describe 
the time evolution of the variational function, the phase 
factor i[3(t)r 2 is needed [9]. The choice of a Gaussian 
shape for the condensate is well justified in the limit of 
weak interatomic coupling, because the exact ground- 
state of the linear Schrodingcr equation with harmonic 
potential is a Gaussian. Moreover, for the description of 
the collective dynamics of Bose-Einstein condensates, it 
has already been shown that the variational technique 
based on Gaussian trial functions leads to reliable results 
even in the large condensate number limit [9]. By in- 
serting the trial wavefunction in Eq. (1), after spatial 
integration, the action of the system reads 
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where the effective potential energy is given by 



W{a) = ^ +3a 2 + 2N(^ + ^-T 29 ( X ^)) , (5) 

with 7 r = (2/7r) 1/2 a r /a , n = (2/7r) 1/2 a (a s - a r )/r 2 , 
x = 2~ 1 / 2 a /r e , r 2 = ag(a s - a r )jr\ and g(x) = 
exp(x)(l — erf(x)), where erf(x) is the error function. 
Note that the action does not depend explicitly on (3, 
which is fully determined by a via the Euler-Lagrangc 
equation (3 = —&/(2<j) (for further details see [10]). 

The extrema of W(a) are obtained as solutions of an 
algebraic equation, which gives the number of bosons as 
a function of the size a of the cloud 
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This equation has either one or three positive roots de- 
pending on the parameters and on number N of atoms in 
the cloud. When three solutions are present, the interme- 
diate one represents an unstable state (i.e. a local maxi- 
mum of the energy), while the other two respectively de- 
scribe a low-density metastable solution (local minimum) 
and a high-density stable solution (absolute minimum) 
within GP approximation. In the local case, Eq. (6) 
reduces to N = (cr 5 — where 7 = (2/7r) 1 / 2 |a s |/ao. 

The variational results for four trap sizes are shown 
in Fig. 1, where the size of the condensate is plotted as 
a function of N (see also Reatto, Parola and Salasnich 
[4]). For comparison, we also show the radius of the cloud 
when a local interaction is assumed. In this case there 
is a critical number N c ~ 0.67ao/a s of bosons beyond 



which the cloud collapses [3]. Fig. 1 shows that the ef- 
fects of non-locality are always important for very tight 
traps while for larger traps non-locality is important just 
when the radius of the cloud rapidly drops for increas- 
ing N. This "transition" is discontinuous for large traps. 
By reducing the trap size, however, this discontinuity is 
strongly reduced and, below about ao = 0.3 fim, the un- 
stable branch disappears and there is a smooth evolution 
from a very dilute cloud to a less dilute state with an 
increasing density as N grows. For large N the size of 
the cloud is independent of the trap size and the atoms 
are in a macroscopic quantum self-trapped (MQST) con- 
figuration. Note that we have also computed the exact 
solution of the GP equation. The variational approach is 
always very close to the exact solution [4] . 

The condensate undergoes density oscillations around 
the minimum of the energy W{a). By performing a 
quadratic Taylor expansion around the minimum one 
finds simple expressions for the monopolc collective fre- 
quency. In the local case we have to — (5 — cr -4 ) 1 / 2 , 
such that lo = 2 for a = 1 (ideal gas), and lo — > for 
(To — » 5 -1 / 4 (i.e. for N — > N c ). Instead, in the nonlocal 
case we obtain 
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where a is related to TV by Eq. (6). For the larger traps, 
where the two branches are present, the frequency of the 
metastable branch is very close to the result given by 
the local approximation. In the stable branch lo starts 
from zero and it raises rapidly as N 1 / 2 . For traps of 
intermediate size, lo has a dip in the transition region 
between the low-density state and the self bound state. 
For very small traps, there is only one branch and the 
frequency increases smoothly with the number of bosons. 
Note that the variational monopole frequency and the 
numerical one (non local case), obtained by solving the 
nonlocal time-dependent GP equation, differ by less than 
3% (see Reatto, Parola and Salasnich [4]). 

We have seen that the low-density and high-density 
branches represent local minima of the effective poten- 
tial energy W{a) of the condensate, given by Eq. (5). 
We know from quantum mechanics that a state that is 
concentrated in one minimum of the potential energy 
may tunnel into a lower one. Ueda and Leggett [11] 
have shown that in the local case, very close to N c , 
there is macroscopic quantum tunneling (MQT) from the 
metastable branch to the collapsed state, which has zero 
radius and energy equal to minus infinity. In our more re- 
alistic picture, there is MQT from the low-density branch 
to the high-density one. The rate of MQT can be calcu- 
lated by the semiclassical formula Tt = Ae~ s / h , where 
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S is the istantonic action, that is the action S with 
imaginary-time (Euclidean action) evaluated along the 
trajectory that makes it extremal [12]. For a quadratic- 
plus-cubic potential the prefactor A is given by A = 
iv(15S B /2-kK) 1/2 [11]. After some tedious but straight- 
forward calculations, one finds 
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where W, N and ui are given by Eq. (5), (6) and (7), 
respectively. In Fig. 2 we plot the rate of MQT for four 
different traps. As expected, this rate becomes large only 
near N c , both in the local and nonlocal case. Moreover, 
Tt increases strongly by reducing the trap size a,Q. As 
previously discussed, in the nonlocal case, below about 
ao = 0.3 /um the unstable branch disappears and there is 
no more MQT. 

We observe that, also if the high density branch is me- 
chanically stable, it has a very short life-time due to 
two-body dipolar collisions and three-body recombina- 
tion. Atoms that inelastically collide acquire substantial 
energy and are ejected from the trap. The total loss 
rate due to two-body dipolar collision and three-body re- 
combination collisions is given by Tl — K J d 3 r|^(r)| 4 + 
L J d 3 r\(j){r)\ 6 , where K = 1.2 • 1(T 2 fim 3 sec" 1 is the 
two-body coefficient and L = 2.6 • 1CP 4 /im 6 sec -1 the 
three-body coefficient [13]. By using the trial wavefunc- 
tion of Eq. (3), we find 
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In this formula the value of TV is related to that of a 
through Eq. (6). In Fig. 3 we plot the total loss rate Tl 
as a function of the number of Li atoms. During MQT 
the condensate density grows rapidly, thereby increasing 
the decay rate from inelastic two- and three-body colli- 
sions. In the case of a unique stable branch, the loss rate 
increases very fast with the number of atoms due to the 
high-density of the system. For sake of completeness, in 
Fig. 4 we show the condensate size, the MQT rate and 
the loss rate for the trap of the Rice experiment (a = 3 
/im). Actually, the trap of the Rice group has a small 
anisotropy and we have used the geometric average of 
the frequencies [1,4]. It is interesting to observe that for 
this trap the high density branch exists also for a number 
of atoms much smaller than N c ~ 1400. 

In conclusion, we have studied the effect of a nonlocal 
inter-atomic interaction for a Bose-Einstein condensate 
of 7 Li atoms in a isotropic harmonic potential at zero 
temperature. We have shown that, in addition to the 
low-density metastable branch, which is present also in 
the case of local interaction, a new stable branch ap- 
pears at higher densities. For a large number of atoms, 
the condensate is in a novel macroscopic quantum self- 
trapped (MQST) state. Very close to the low-density 



critical threshold, there is macroscopic quantum tunnel- 
ing (MQT) between this low-density metastable branch 
and the high-density one. This prediction is more realis- 
tic than that of MQT from the metastable branch to a 
collapsed state, which has zero radius and energy equal 
to minus infinity. We have calculated the rate of MQT 
with the semiclassical technique of istantons. Finally, we 
have considered the decay rate of the Bose condensate 
due to inelastic two-body dipolar collisions and three- 
body recombination. The life-time of the condensate in 
the high-density branch is very short because of the small 
radius of the condensate in these configurations. Due to 
the very short life-time of 7 Li, experimental studies of the 
high-density phase of the condensate are practically im- 
possible. Nevertheless, this is only a technical difficulty, 
perhaps surmountable by using a different atom, and it 
is imaginable that in the future our predictions can be 
checked. 
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FIG. 1. Size of the condensate as a function of the number 
N of 7 Li atoms. Four different traps, from top to bottom and 
from left to right: ao = 1.0 /im, 0.6 /im, 0.3 /xm, 0.1 /im, 
respectively. The dashed lines represent the results with local 
interaction. 
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FIG. 2. Rate of MQT as a function of the number N of 7 Li 
atoms. Four different traps, from top to bottom and from left 
to right: ao = 1.0 /im, 0.6 fim, 0.3 /an, 0.1 /im, respectively. 
The dashed lines represent the results with local interaction. 
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FIG. 3. Loss rate due to two- and three-body collisions as a 
function of the number N of 7 Li atoms. Four different traps, 
from top to bottom and from left to right: ao = 1.0 fim, 0.6 
/tm, 0.3 fim, 0.1 /tm, respectively. The dashed lines represent 
the results with local interaction. 



4 



3.0 



E ' 
3, 

b 1.0 - 



<D 

W 400 - 
W 

E 
o 

S 200 - 



E ' u 
o 

"1 10° 



FIG. 4. Condensate size a, MQT rate Ft and loss rate Tl 
as a function of the number N of 7 Li atoms. Trap of Rice 
University with ao = 3.0 /jm [1]. The dashed lines represent 
the results with local interaction. 
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